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Abstract. We prove the existence of a global bifurcation branch of 
27r-periodic, smooth, traveling-wave solutions of the Whitham equation. 
It is shown that any subset of solutions in the global branch contains 
a sequence which converges uniformly to some solution of Holder class 
C a , a < i. Bifurcation formulas are given, as well as some properties 
along the global bifurcation branch. In addition, a spectral scheme for 
computing approximations to those waves is put forward, and several 
numerical results along the global bifurcation branch are presented, in- 
cluding the presence of a turning point and a 'highest', cusped wave. 
Both analytic and numerical results are compared to traveling-wave so- 
lutions of the KdV equation. 

1. Introduction 

The Whitham equation, 

Vt + ^VVx + K ho *r) x = 0, (1.1) 

combines a generic nonlinear quadratic term with the exact linear dispersion 
relation for surface water waves on finite depth. Here, the kernel Kh ■ = 
l (ch ) is the inverse Fourier transform of the phase speed 



c , (e) : =y^p2i (L2) 

for the linearized water-wave problem; the constants g, ho and en := y/gho 
denote, respectively, the gravitational constant of acceleration, the undis- 
turbed water depth, and the limiting long-wave speed. The function rj(t, x) 
describes the deflection of the fluid surface from the rest position at a point 
x at time t |19] , The Whitham equation was introduced by Whitham (see 
|18] ) as an alternative to the Korteweg-de Vries (KdV) equation, which de- 
scribes the evolution of predominantly uni-directional, small-amplitude sur- 
face waves in shallow water [21 El \TU\ . The KdV equation appears from (|1.1|) 
if one replaces c/j (£) with its second-order approximation at £ = 0, 

co - \c Q hlf, (1.3) 

and this explains the main motivation behind the Whitham equation. In 
particular, since £ = x may be interpreted as a wave number, A being a 
typical wavelength, one sees from (jl.3p that the KdV equation is a poor 
approximation to (|1.2p for large wave numbers. 

The Whitham equation (jl.ip with the kernel (|1.2p has some very inter- 
esting mathematical features. First of all, it is generically nonlocal, mak- 
ing pointwise estimates difficult. Moreover, c/ lQ (£) has slow decay, and the 
kernel is singular (it blows up at x = 0). This makes the Whitham 
equation in some important respects different from many other equations of 

l 
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the form (jl.ip . One example is the possibility of singularities in solutions: 
the Whitham equation features (phenomenological) wave-breaking |15j 
i.e. bounded solutions with unbounded derivatives — and is conjectured to 
admit solutions with cusps |18| . 

In this paper we consider steady solutions of the Whitham equation, i.e., 
traveling-wave solutions characterized by a constant speed and shape. The 
existence of smooth, small-amplitude, periodic traveling-wave solutions was 
established in |10j . and their properties numerically investigated. The nu- 
merical calculations pursued in that investigation indicate a global branch of 
solutions approaching a 'highest', cusped, wave — as well as a sequence of 
periodic waves converging to a solitary wave. An analytic proof of the latter 
fact is given in [9], for small waves, using a variational approach developed 
mainly in [12] . The aim of this note is: (i) to place the Whitham equation 
in a general functional-analytic framework, in which the bifurcation theory 
for that and many related equations is quite natural; (ii) to extend the main 
local bifurcation branch found in |10| to a global one; and (iii) to provide a 
spectral scheme for efficient calculation of Whitham waves. In addition, we 
prove some convergence results along the global branch of solutions, and sev- 
eral numerical calculations are given, complementing the more abstract exact 
theory. In particular, numerical evidence is given for a turning point along 
the global bifurcation branch, as well as additional calculations supporting 
the existence of a cusped wave at the end of the bifurcation branch. 

Judging from the forms of the KdV and Whitham equations, one would 
expect several similarities for small waves of small wavelength (one such con- 
vergence result is given in [9]). To illustrate this, we perform a comparative 
local bifurcation analysis for the two equations, analytically and numeri- 
cally. Although the solutions arise in the same manner from a curve of 
vanishing solutions, they develop differently along the corresponding bifur- 
cation branches. We show that the Whitham solutions are all smooth and 
subcritical (their normalized wave speed is uniformly bounded away from 
unity) , and that they converge uniformly to a wave of C Q -regularity, a < | . 
The intricate form of the Whitham kernel has yet hindered us from proving 
that the limiting wave is cusped and non-trivial. To see the difficulty, note 
that the Whitham kernel Kh is not easily seen to be positive definite (cf., 
e.g., HU). 

This is the disposition of the paper: in Section [2] we normalize the equa- 
tion for traveling wave-solutions and briefly discuss Fourier multipliers on 
Holder spaces. Section [3] contains a functional- analytic formulation of the 
problem, and local bifurcation theory for the Whitham and the KdV equa- 
tions. Section [J] is the main analytic section, in which we state and prove the 
global results. The sections [5] and [6] are devoted to the numerical analysis. In 
Section [5] we describe a spectral scheme for the Whitham equation, and use 
it to calculate wave profiles. Finally, Section [S] takes a look at the shape of 
the bifurcation branches, and provides a numerical comparison of the KdV 
and Whitham waves. 
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2. Some preliminaries 

Throughout the paper, the operator T will denote the extension to the 
space of tempered distributions S'(M) of the Fourier transform 

Hf)(0 ■= J f(x)exp(-ixOdx 

on the Schwartz space 5(M), with inverse J-~ 1 (f)(x) := ^ f /(£) ex.p(ix£) d£. 

TRAVELING WAVES. By considering steady solutions with wave speed c > 0, 

tp{x — ct) := r](t, x), 

the left-hand side of (jl.ip may be integrated to — op + ff^</? 2 + Kh *<fi- bi the 
following, we keep with the convention in denoting the steady variable x — ct 
again by x. The scalings ■^ l f > <p and x i— > h^x then yield the normalized 
equation 

- fi<p + ip 2 + K*<p> = 0, (2.1) 
where /i := c/cq is the non-dimensional wave speed, and we have set the 
constant of integration equal to zero (see Remark 12.11 below for a discussion 
of this) . The symbol of the convolution operator in (|2.ip is given by 

F{K)(t):=(^) 1 ' 2 (2.2) 

for the Whitham equation. The corresponding Fourier multiplier for the 
KdV equation has the form 1 — g-£ 2 , and the equation reads 

(1 - flip + ip 2 + \<p>" = 0. (2.3) 
o 

Any solution of (|2.ip with the kernel given by (|2,2p will be called a Whitham 
solution, and any solution of (|2.3p will be called a KdV solution. We shall 
be dealing only with 2-7r-periodic solutions. 

REMARK 2.1. The normalization B = in 

— [up + tp 2 + K*ip = B, 

cf. (|2.ip . is a choice of convenience: both the Whitham and KdV equations 
are invariant under the Galilean transformation 

<p i-)- (p + 7, /u i->- fi + 27, B i->- B + 7(1 — /j, — 7), 

for any 76R. As an alternative to letting B = one may therefore consider 
waves of vanishing mean amplitude. 

Fourier multipliers ON Holder spaces. The analysis given in this pa- 
per relies on certain properties of Fourier multiplier operators given by clas- 
sical symbols, and we briefly recall some of those. A smooth, real-valued 
function a on R is said to be in the symbol class S m if for some constant 
c > and any non- negative integer k, the estimate 

\dla(0\<c(l + \£\) m - k 

holds. Analyzing the action of such operators on the Sobolev spaces H s is 
straightforward owing to Parseval's formula. However, the analysis is some- 
what more subtle if classes of continuously differentiable functions or other 
Banach spaces are considered. Although there exist natural scales of function 
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spaces for this class of operators (cf. |17j ) we have chosen to work with the 
classical spaces C| ven (§), s > 0, consisting of 27r-periodic even functions that 
are times continuously differentiable with the [sj th derivative bounded 
and continuous, and Holder continuous with Holder exponent s — [s\ for 
s Z. We then have the following lemma. 

LEMMA 2.2 ([]]). Fix m 6t, a e S m , and let s <E M+ \ Z be a non-integral 
non-negative real number. Then 

a(D) G £(C(S) s+m ,C7(S) s ) 

is a bounded linear operator of order —m. 

REMARK 2.3. If one extends instead the family of Holder spaces C s , for 
s Z, with the so-called Zygmund spaces for integral values of s, one attains 
a scale of spaces for which Lemma [2 .21 is valid also for integral values of s (for 
more information on the relations between these different function spaces, 
see |17|). The introduction of Zygmund spaces is, however, not necessary to 
obtain the results in this paper. 

3. Comparative local bifurcation theory 

From now on, let ^ < a < 1. We consider C^ en (S), the space of even and 
a-H61der continuous real- valued functions on the unit circle §. The Whitham 
equation contains a generic nonlocal smoothing operator in the form of the 
Fourier multiplier (tanh(^)/^) 1 / 2 , whereas the corresponding KdV multiplier 
1 — g£ 2 is a second order differential operator. This difference may be easily 
overcome: rewriting the KdV equation in non-local, smoothing, form enables 
a similar treatment of the two equations. To illustrate how the analysis used 
for the Whitham equation can be applied to a larger class of equations, a 
comparative local bifurcation analysis is performed for the Whitham and 
the KdV equation (the KdV equation being a cardinal example of dispersive 
evolution equations and the Whitham equation's closest kin). 

THEOREM 3.1 (Functional-analytic formulation). Fix a > ^ and let fi > 
0. The solutions in C^ en (§) of the Whitham equation (|2.ip . and those of 
the KdV equation (|2,3p . coincide with the kernel of an analytic operator 
F: C c Q vcn (S) x R >0 -> C c Q vcn (S) given by 

F(tp, n) := fop - L(fj,)tp + N(ip, fi), 

where L(^i) is bounded, linear and compact, and 

Nw(tp, M) = -^ 2 > N Kdv (ip, fi) = - L(fi)tp 2 , 

fulfills D^N[0,^] = 0. Thus D^F[0, /i] is Fredholm of index 0. In the case 
of the Whitham equation the operators L and N are independent of fi. 

Proof. We put the KdV equation in the form 

where (l — r 2 5 2 ) _1 / := exp y — ipj*f is the inverse of l-r 2 <9 2 on <S'(R), 
the Schwartz space of tempered distributions on 1R. Define 

L w := K * and L Kd v := (l - ^<9 2 ) . (3.1) 
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For functions / £ C^ ven (§) with a > 1/2 one has 

f(x) = a/; cos(A;x) and [g^l < oo, 

fe>o fc>o 

meaning that C^ en ) is a subalgebra of the Wiener algebra of 27r-periodic 
functions with absolutely converging Fourier series [13]. For such functions 
rudimentary calculations reduce (|3.ip to the intuitive formulas 

L w f(x) := £ a k (^^) V2 cos(fex), (3.2) 
fc>0 

L KdV (Li)f(x) := Y : ° i , 2 cos(fcr), (3.3) 

where for the Whitham equation one must use the fact that K £ L X (R) (for 
details on if and its symbol, see |10j). The Fourier multiplier symbols in (|3.2p 
and ([33]) belong to the symbol classes S~ 1/2 (M) and S"~ 2 (IR), and §Z/2§ and 
(|3.3p are therefore bounded linear operators C" ven (§>) — > C^^(S) for s = 1/2 
and s = 2, respectively, see Lemma [2.21 Each of them is also invertible with 
bounded linear inverse L _1 : C^T*(S) — > C^ en (S) for the same values of s. 
Due to the compactness of the embedding Ce ve nC§>) C^ ven (§i), /3 > a, both 
operators are compact on C" ven (§>). We may thus define mappings F w and 
F KdV : C« ven (S) x M >0 ->• Cg en (S) by 

F w ((/5,/i) := /x^ - L w v? - (^ 2 , (3.4) 

F K dv(^, M) := - L Kdv (/i)<p - L Kc iv(m)v 2 - ( 3 -5) 

In both cases .F is analytic in both its arguments (L KdV (//) is the inverse 
of 1 — ^<9 2 , which is analytic in /i). We also have F(0,li) = 0, and the 
linearization F[0, fi] = /iid — L is Fredholm of index (it is a continuous 
perturbation of the identity by a compact operator). □ 

Local bifurcation. 

COROLLARY 3.2 (Local bifurcation, Whitham). 

(i) SUBCRITICAL BIFURCATION. For each integer k > 1, there exist fi k : = 
(tanh(/c)//c) 1 ^ 2 and a local, analytic curve 

mM e ),^))eCcn(S)x(0,l) 

of nontrivial 2n/k -periodic Whitham solutions with D E f(0) = cos(kx) that 
bifurcates from the trivial solution curve fi >— > (0, /u) at (<^(0), fJ.(0)) = (0,/Xfc). 
In a neighborhood of the bifurcation point (0,/ifc) these are all nontrivial solu- 
tions ofF\y((p,fi) = in C^ ven (S) x (0, 1), and there are no other bifurcation 
points fi > 0, n 7^ 1, for solutions in C^ ven (S). 

(il) TRANSCRITICAL BIFURCATION. At fi = 1 the trivial solution curve 
H i — V (0,/u) intersects the curve /i 4 (fj, — 1, /i) o/ constant solutions ip = 
fx — 1; together these constitute all solutions in C^ cn (S) in a neighborhood of 
(<p,H) = (0,1). 
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Proof. The fact that F w [0, fi) = /iid — L w is Fredholm of index and the 
formula (|3.2p show that are all simple eigenvalues of L w , and that no other 
eigenvalues fj, > exist. The assertion then follows from the analytic version 
of the Crandall-Rabinowitz theorem for bifurcation from a simple eigenvalue 
[6l Thm 8.4.1]. The fact that the solutions are 27r//c-periodic can be seen by 
restricting attention to the subspaces {(p £ C" ven (S): (f is 27r/A;-periodic}, 
and, corresponding to the case k = 0, it is instantly verified that <p = fj, — 1 
is a solution. By uniqueness, this family must therefore constitute the local 
bifurcation curve at fi = 1. Since for all other [i > the linearization 
F w [0, n] is Fredholm of index zero with trivial kernel, it is a consequence 
of the implicit function theorem that the vanishing solution is, locally, the 
unique solution in C" ven (§>). □ 

REMARK 3.3. A more detailed proof of local bifurcation in a somewhat larger 
space, for k = 1, but with arbitrary wavelength, was given in |10| . 



For the KdV equation the change of variables 

ip^k 2 ip(k-), ^^fc 2 (/i-l) + l, (3.6) 

means that the local bifurcation analysis can be reduced to the case k = 1. 

COROLLARY 3.4 (Local bifurcation, KdV). Let fi* := 5/6. A local, analytic 
curve e t- > (99(e), //(e)) E Cg ven (S) x M>o of nontrivial KdV solutions with 
D £ (f(0) = cos(x) bifurcates from the trivial solution curve at (0, fj,*). 

Proof. We have that 

L\,F K dv[0,/i]v? = ^a fc (fi - -^-p ) cos(kx). 

Therefore, ker (D^ F K dv[0, /i*]) = span(cos(x)) and 

D ^ F Kdv[0, (J.*] cos(x) = § cos(x) ran (D^, F Kd v[0, M*D • 

As in the proof of Corollary 13. 2| since F K dv[0, fJ-*] is Fredholm of index 
0, the assumptions of the Crandall-Rabinowitz theorem are fulfilled. □ 

REMARK 3.5. The bifurcation points corresponding to k = 1 are not the 
same for the Whitham and the KdV equation: we have fi\ ~ 0.87 but 
H* » 0.83. 



REMARK 3.6. In the original variables, before the transformation (|3.6p . the 
wave speed \i of the KdV equation becomes negative for k > 3, a fact related 
to that KdV is a proper model for small water waves when the wave speed 
is positive and the wave number small. 

4. Global bifurcation for the Whitham equation 

As earlier, let a > \ denote a fixed Holder exponent. Let F w be the 
Whitham operator from Theorem 13.11 defined by (|3.2p and (|3.4p . With 

U := {(if, fj) e C e Q ven (§) X (0, 1): if < , 

we let 

S-j^^Gt/: F w ( V ,/ t )=0} (4.1) 
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be our set of solutions. By C(X) we denote the Banach algebra of bounded 
linear operators on a Banach space X, and by Ci s (X) those which also have 
an inverse in C(X), i.e., those which are Banach space isomorphisms. 

LEMMA 4.1 (L°°-bound). Let fi > 0. Any bounded Whitham solution satisfies 

IMloo < M + || Lw|| j c(l«(s))- ( 4 - 2 ) 
Proof. From flip — L w — ip 2 = we obtain that 

M 2 < mM + II W ||£(l°°(s))IMIoo, 

where we have used the fact that the Whitham kernel is integrable. Either 
if = 0, or we may take the supremum and divide by ||y||oo- in either case, 
(IP) holds. □ 

LEMMA 4.2 (Fredholm). The Frechet derivative D^FwIv^/i] is a Fredholm 
operator of index for all (ip, fi) G U. 

Proof. We have 

D^, F w [tp, (j] = (fjL- 2tp) id - L w , 
and, for any given (cp, [/,) G U, that (/i — 2ip) id G £i s (C^ en (S)). In view of 
that L w is compact on C Q (S), the operator D^, F w [</?, /i] is Fredholm. The 
linearization D (/5 F w [0, f/] has Fredholm index zero along the trivial solution 
curve; we have 

T^(fM-2Tip) id-Lw G<7([0, l],£(C a (§))), 

and since the index is continuous in the operator-norm topology, it follows 
that it is zero also at (ip,fi). □ 

LEMMA 4.3. Whenever ((p, fi) G S the function ip is smooth, and bounded 
and closed sets of S are compact in C^ ven (S) x (0, 1). 

Proof. We write the Whitham equation in the form 

^ = F(^, / .):=|-(^-L w ^ 1/2 . (4.3) 

The mapping L w is bounded and linear C a (§) -> C a+1 / 2 (S,) (cf. Eemma lz.zl ). 
and x i — y \fx is real analytic for x > 0. Consequently, if we let 

V:= {(<p,fi) GC Q (S) x (0,1): ^ > L w 9?} , 

then F is real analytic V -> C Q+1 / 2 (§). The space C Q+1 / 2 (§) is relatively 
compact in C Q (S), whence F maps bounded subsets of V into pre-compact 
sets. We may then prove: 

Smoothness. For any ip G S there exists a constant R\ such that sup y> < 
R\ < fi/2. Since ip is a fixed point of F(-,/i) we have (<p,f£) G V. A straight- 
forward induction argument reveals that ip G C°°(S). 

COMPACTNESS. Let if C S 1 be bounded and closed in the C a (§>) x R- 
topology. Then K C V, and {93: (ip,fx) G if} = Fif is pre-compact in 
C a (§). Any sequence {(<^j, Mj)}j>i C K thus converges to a pair (</?o,/io) i Q 
the C a (S) x R-topology. The fact that K is closed implies that (</?o, Mo) G K, 
whence K is compact. 

□ 
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Global BIFURCATION. Using Lemmata [O] and [O] we shall now show that 
the local branches of Whitham solutions in Corollary 13.21 can be globally 
extended. The following result is immediate (cf. [6]) if we are able to show 
that, in some small neighborhood |e| < 5, //(e) is not identically equal to a 
constant. 

Theorem 4.4 (Global bifurcation). The local bifurcation curves e t— > ((p(e), fi(e)) 
of solutions to the Whitham equation from Corollaru VS. OA extend to global con- 
tinuous curves of solutions R>o — S- S, with S as in (|4.ip . One of the following 
alternatives holds: 

(i) || y(e) || era (g) is unbounded as e —> oo. 

(ii) The pair (y(e),/t(e)) approaches the boundary of S as e — >■ oo. 

(iii) The function e *— > (</?(e), /t(e)) is T -periodic, for some T £ (0, oo). 

Proof. According to [6l Thm 9.1.1], and in view of Lemmata 14.21 and 14. 3| the 
assertion follows if any of the derivatives n^ k \0) ^ 0. It will be shown in 
Theorem |0] that /i(0) / 0. □ 

In order to establish the bifurcation formulas for Whitham equation — and 
thereby Theorem 14.41 — we apply the Lyapunov-Schmidt reduction (cf. |14j). 
Let /t* := /ii be the bifurcation point from Corollary 13.21 and let 

(f*(x) := cos(x). 

Let furthermore 

M := i ^a fc cos(A:x) G C7 a (§) I , 

W 1 J 

and 

N := ker (D^ F w [0, fi*]) = spanfa*). 

Then C^ en (S) = M © N and we can use the canonical embedding C a (§) 
L 2 (S) to define a continuous projection 

with (u,v) L 2 {s) : = i f § uvdx. 

THEOREM 4.5 (Lyapunov-Schmidt, |14j). There exists a neighborhood O x 
Y CU around (0,//*) in which the problem 

F w { ( p,n) = (4.4) 

is equivalent to that 

$(e^*,/z) :=nF ff (^ + ^,^/i) = (4.5) 

for functions ip £ C°°{O n x Y,M) ; $ £ C°°{p N x YjJV), and Ojy C JV 
an open neighborhood of the zero function in N. One has $(0,/i*) = 0, 
V>(0,//*) = 0, Dm "0(0, /t*) = 0, and solving the finite- dimensional prob- 
lem (|4.5p provides a solution ip = eip* +ip(eif*, fi) of the infinite- dimensional 
problem (|4.4p . 
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Theorem 4.6 (Bifurcation formulas). Let fi* = ^/tanh (1), C\ := (/i*-l) _1 , 
and C 2 := (2/i* - y/2 tanh (2)) _1 . T/te mam bifurcation curve (k = 1) for 
the Whitham equation found in Corollary \3.S\ satisfies 

ip(e) =ecos(x)+e 2 (±Ci + C 2 cos(2x)) +C(e 3 ), (4.6) 

and 

M ( e ) = fl * + £ 2 (C l + C 2 ) + 0(s 3 ), (4.7) 

m C^ en (S) x (0,1) as e — > 0. In particular, /i(0) < 0, Corollary VS.'A (i) 
describes a subcritical pitchfork bifurcation, and Theorem\4-4\ holds. 



REMARK 4.7. The solution curve <p(e) for the Whitham equation differs 
significantly from that of the KdV equation. Indeed, relying on Corollary 13. 41 
it is possible to calculate the asymptotic expansion 

^Kdv(e) = e cos(x) + e 2 (cos(2x) - 3) + C(e 3 ), 

in C" ven (S) as e — )• for the periodic KdV-solutions found there. Although 
this agrees with (j4.6|) to the first order in e, it does not for e 2 . In particular, 
the coefficient in front of cos(2x) is positive for the KdV-solutions, whereas 
it is negative for the Whitham solutions. 

Proof. We perform first the analysis for \i. We already know that e i— > //(e) 
is analytic at e = and that /i(0) = //*, so it remains to show that /i(0) = 
and to determine (1(0). We have 



D^F w [0,/i*](^,^) = -2^ 2 



D^F w [0,^]^ = ^, 

and the value of (i(0) may be explicitly calculated — for the bifurcation for- 
mulas used in this proof, we refer to |14l Section 1.6] — as 

. ( ) = l <D^F w [0,/x*](^,(p*),<p*) L 2( S) = 
2 (D 2 ^F w [0,^]ip*,^) L 2 {s) 

since Jg cos 3 (x) cfe = 0. 

Moreover, when fi(0) = one has that 



1 



A(0) = /r>2 = [n , „, ^- (4.8) 

Since D 2 ^ F w [0, //*] = id we find that the denominator is of unit size. One 
calculates 

= F w [99 + ip(ip, n),fi](<p* + D v ip[ip, /j]ip*,(p* + D^, /j]ip*) 

+ UB ip F w [Lp + i/;(<p, fj),n] B 2 W i/j[<p, n](<p*,<p*), 
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and, in view of that F w is quadratic in <p, 

= 3 n D ^ F w [ip + t%, n) , n\ ((p* + D v i> [if, iA <p* , D J„ il>[<p, fi] ,<p*)) 
+ n D v Fw[y + /*),/•*] </#, n\ (<p* ,<p*,<p*). 

Using the form of D^Fw together with that tp(0,fi*) = ip[0, fi*]ip* = 
one finds that 

D J w *[0, n (m* id - Lw) w ^[0, /x*] , <p* , ^*) 

-en>* D^v[o,/i*](v?*,^). 

We have ran(/x* id — L w ) = M, so that IT (/x* id — L w ) = 0. We thus need 
to determine ip* ^[0, /J,*](tp*, if*). Since D^FwfO, /x*] = \i* id — L w is an 
isomorphism on M, it is possible (see again |14| Section 1.6]) to rewrite 
D^[0,xx*](^*) as 

D^^[0,/x*](^,^) 

= - (D v F w [0, /x*])" 1 (id - n) D w F w [0, At*](v>*,¥>*) 
= - (D„ F w [0, /x*])" 1 (id - II) ( - 2<p* 2 ) 
= (D^F w [0,/i*])~ 1 (l + cos(2x)) 
cos(2:e) 



A** - 1 /x* - v / tanh(2)/2 
After multiplication with cos(x) this equals 

cos(x) cos(x) cos(3x) 



A** - 1 2(/i* - v / tanh(2)/2) 2(/x* - v / tanh(2)/2) 

In view of (|4.8|) and (|4.9p the coefficient in front of cosx equals |/x(0). All 
taken into consideration, we obtain (|4.7p via a Maclaurin series, and one 
easily checks that /x(0) < 0. 

To prove (|4.6p one makes use of the formula 

<p(e)=E<p* + il>{e<p*, l i.{e)) (4.11) 

from the Lyapunov-Schmidt reduction (cf . Theorem 14. 5p . We already know 
that <p(0) = and <p(0) = cos(x), so it remains to calculate </5(0). It follows 
from fltTTU) that 

(p(e) = B 2 W m n*](<p*,<p*) + 2 V>[0, xx *](</, £(0)) 
+ V[0, /x*](/i(0), A(0)) + ^[0, xx*]xx(0). 

Since -0(0, /i) = where -i/> exists, we have D„^(0, xx*) = 0. Combining this 
with fi(0) = one finds that 

= D^^[0,jLi*](cos(x),cos(a;)), 

so that the proposition now follows from (|4.10p . □ 
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Properties along the bifurcation branch. 

THEOREM 4.8 (Uniform convergence). Any sequence of Whitham solutions 
(tprnf^n) £ $ h as a subsequence which converges uniformly to a solution 
(f £ C a (§>), with a £ (0, ^) arbitrary. At any point where cp < ^ the solution 
is a-Holder continuous for all a £ (0, 1). 

Proof. Since \i n £ (0, 1), Lemma [4. 1 1 implies that ip n is uniformly bounded in 
C(S). The proof of Theorem 4.1 in |10] shows that any uniformly bounded 
sequence of Whitham solutions (i.e. in L 00 (R)) is equicontinuous; this follows 
from a standard argument using the integral kernel of L w . Since we are 
dealing with periodic solutions, we may apply the Arzola-Ascoli lemma to 
conclude that a subsequence of ip n converges in C(S). 

Since L w maps C(S) into C 1//2 (§) (cf. Lemma [2 .2p we see from (|4.3p that 
(p € C 1 / 2 (S) wherever 2(p ^ fi. On the other hand, when ip(x) = ^ we have 
that 

A* /V V /2 

\ip(x) -tp{y)\ = --(f(y) = f — -L w <p(y)\ 

= (L w <p{x) - L w <f(y)) 1/2 < C\x - y\ l l\ 

in view of that L w V 9 i s a-Holder continuous for a=\. Hence ip £ C 1//4 (S). 
Repeating this argument once shows that ip £ C 3 / 4 wherever 2(p ^ fi, and 
ip £ C 3 / 8 (§). Bootstrapping yields Holder continuity below a = 1 for 2y> ^ 
/i, and below a = | for tp £ C a (S). □ 



PROPOSITION 4.9 (Characterization of blow-up). Alternative (i) in The 
rem \4-4\ can happen only if 



liminf inf ( ^ - <p(x;e) ) = 0. (4.12) 



e— >oo xd 



In particular, alternative (i) implies alternative (ii). 

Proof. Assume that lim inf £ _>.oo inf x — V 7 ( a; ) e )^ — ^ f° r some 5 > 0. 

Any such solution of the Whitham equation satisfies 

i / \ / \ I \L w tp(x) -L w ip(y)\ | L w ip(x) - L w <p(y)\ 

\(p{x) - cp{y) = j- —— < — . 

H - (p(x) - ip{y) 25 

Since L w is continuous C(S) — > C 1 / 2 (S) and the family {(p(e)} £ is uni- 
formly bounded in C(S) (cf. LemmaUTJ, it follows that {(p(e)} £ is uniformly 
bounded in C 1 / 2 (S) too. Repeating the argument for L w as a continuous 
operator C 1 / 2 (S) -> C a (S), a < 1, yields that 

llv(e)llo-(S) <C<T 2 , «e(0,l), 

for some constant C depending only L w (recall that fi is bounded by as- 
sumption). This shows that || 99(e) ||c a (S) — ?• 00 is possible only if (|4.12p 
holds. □ 
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SUMMARY. In Section U] we have shown that each of the subcritical bifur- 
cation branches in Corollary 13,21 (i) can be extended to a global bifurcation 
branch in the set S of solution pairs (</?, /x) G C" ven (§) x (0,1) satisfying 
(f < n/2. These solutions are all smooth. Each bifurcation curve contains a 
sequence of solutions which converge in C° to a solution of Holder regularity 
a £ (0, i), with the better pointwise regularity a £ (0, 1) at any point x 
where ip(x) < fi/2. If the C Q -norm, a > |, blows up along the bifurca- 
tion branch, the corresponding sequence of solutions approaches ip = ^ in 
the sense of (|4.12p . In spite of this — and in spite of the numerical evidence 
available — an analytical proof of a 'highest', cusped, limiting wave at the 
end of each bifurcation branch remains an open problem. 



5. Spectral schemes for the Whitham equation 

In this section, the numerical approximation of solutions of (|l.ip is in view. 
Both for the approximation of traveling waves and for the discretization of 
the evolution problem, spectral schemes are used. This is a natural choice 
since the operator defining the linear part of the equations is given in terms 
of Fourier multipliers. 



TRAVELING WAVES. Solutions of (|2,ip are approximated using a spectral 
cosine collocation method. To define the cosine-collocation projection, first 
define the subspace 



Sn = span^ |cos(/x) 



< / < N - 1 



of L 2 (0,7r), and the collocation points x n = K 2n N 1 for n = 1,...,N. The 
discretization is defined by seeking (j)^ in Sn satisfying the equation 

-H<t> N + ^ N + K N <t>N = ^ (5-1) 

where the operator K N the discrete form of K defined in (|2.2p . The discrete 
cosine representation of 4>n given by 

N-l 

4>n(x) = w (1)®n{1) cos(Zx), 

1=0 

where 

and $at(/) are the discrete cosine coefficients, given by 

N-l 

$at(Z) = w(l) 4>N{xn)cos(lx n ), for I = 0, ...,N - I. 

n=0 

Now if the equation (|5.ip is enforced at the collocation points x n , the term 
K N (pN may be practically evaluated with the help of the matrix [K N ] (m, n) 

by 

N 

[K N ] 4> N (x m ) = [K N ] {m,n)(j> N {x n ), 

n=l 
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FIGURE 1. An approximate traveling-wave solution 4>m with 
nondimensional wavespeed /i = 0.789, wavelength 2-7T, and nondi- 
mensional waveheight 0.3368 is shown in the left panel. The col- 
location approximation used 16 gridpoints, which are indicated as 
dots in the graph. The amplitudes of the discrete Fourier modes 
are shown on the right. 

where the matrix [_K" ] (m, n) is defined by 



[K N ] (m,n) = u> 2 (fc)yj tanhJ cos(7x n ) cos(7x m ). 

1=0 

Thus, equation (|5.1|) enforced at the collocation points x n yields a system 
of N nonlinear equations, which can be efficiently solved using a Newton 
method. The cosine expansion effectively removes the singularities of the 
Jacobian matrix due to translational invariance and symmetry of the solu- 
tions. The nondimensional speed [i is used as the bifurcation parameter, 
but as the turning point is approached, it is more convenient to change to a 
waveheight parametrization. The computation can be started for \x close to 
but smaller than the critical speed (if-, and with an initial guess for the first 
computation given by comparing (|4.6p and (|4.7p . 

In Figure [TJ a periodic traveling wave with speed fi = 0.789 is shown. In 
panel (a), the function is drawn on the interval [0, 2ir], and the collocation 
values are indicated with dots in the curve. In panel (b), the values of the 
discrete cosine coefficients are shown. This computation was done with only 
16 modes. In Figure [2l a periodic traveling wave with speed \jl = 0.7715 is 
shown. It appears that the wave is steeper, and higher cosine coefficients are 
important. The computation was done with 64 modes. To make sure that 
the computed functions are approximate traveling waves for the Whitham 
equation, we have also used a dynamic integrator for the time-dependent 
Whitham equation, as shown in the next subsection. 

TIME-DEPENDENT DISCRETIZATION. For the purpose of approximating time- 
dependent periodic solutions of (jl.ip . a Fourier method is optimal. To define 
the Fourier-collocation projection, the subspace 



•Sn = spanc < exp(ifcr) 



k£Z, -N/2 <k< N/2 - 1 \ (5.2) 
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FIGURE 2. An approximate traveling- wave solution 4>m with 
nondimensional wavespeed fi = 0.7715, wavelength 2-7T, and nondi- 
mensional waveheight 0.4152. The collocation approximation used 
64 gridpoints, and the amplitudes of the discrete Fourier modes 
are shown on the right. 

of L 2 (0, 2ir) is used. The collocation points are denned to be Xj = ^jf- for 
j = 0, 1, TV — 1. Note that N here is in general different from the iV used 
in the previous section. Let P/v be the projection operator onto Sn, and let 
In be the interpolation operator from C^ r ([0, 2ir]) onto 5/\r. As explained in 
[31 IT], this operator is defined for a given u E C^ r ([0, 2ir]), by letting I^u be 
the unique element of Sm that coincides with u at the collocation points Xj. 
The same scaling as in the steady case is used, so that the non-dimensional 
form of the equation is 

r] t + 2r]r] x + K * T] x = 0. (5.3) 

The spatial discretization is then defined by the following problem. Find a 
function r]N '■ [0,T] — > Sn, such that 

/ d t r] N + d x I N {if N ) +1C N * d x r) N = 0, i£ [0, 2vr], , , 

\ w(-,0) =PnVo, { ' 

where 770 is the initial data. If r]j\[(-,t) is written in terms of its discrete 
Fourier coefficients fj^ik) as 

Vn(x) = fj N (k)ex.p(ikx), 

-N/2<k<N/2-l 
the operator K can be evaluated using the formula 

fc N Vn(x) = ^at(O) + ^2 \J i tanh k rj]y(k) exp(ikx). 

l~N/2<k<N/2~l 

Thus the operator fC N is the truncation at the JV/2-st Fourier mode of the op- 
erator given by the periodic convolution with K. Note that this formulation 
includes the truncation of the Fourier mode (Pn(—N/2) which ensures that 
the solution stays real, and which otherwise can lead to instabilities in the 
computation. The time integration of fj5.3j) may be carried out in a number 
of ways. We have chosen to use a midpoint method for both the linear and 
the nonlinear term. This method can be motivated as follows. Suppose the 
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TABLE 1 . Error in evolution code after 5 and 50 periods for the 
traveling wave shown in Figure [T] 



N 


St 


L^-error 


I u\ oo 


— \un\od 


L^-error 


\u\oo — \un\oo 


2 b 


2 -lU 


2.675e-04 


2i 


538e-05 


9.635e-04 


2.612e-05 


2 6 


2-14 


2.644e-05 


li 


\76e-06 


5.299e-05 


1.587e-06 


2 7 


2 -19 


1.021e-06 


5i 


576e-08 


4.489e-07 


5.286e-08 








1 period 




100 


periods 



discrete solution is known at the time t. Then by the fundamental theorem 
of calculus, it is clear that 

/t+St 
d t r](xj,s)ds. 

Now using the equation (|5.4p and approximating the integral with help of 
the trapezoidal rule yields 

rj N (xj,t + St) = r] N (xj,t) + y {d x I N {rf N (xj,t) + rf N (xj,t + St))} 

+ j{lC N VN (x 3 ,t) + IC N rj N ( Xj ,t + 5t)} + (St 3 ). 

From this equation, one may easily construct a fully discrete approximation. 
This can be done either in physical space, or in Fourier space. Constructing 
the approximation in Fourier space appears to be more convenient for the 
time- dependent problem. Since nonlinear terms are involved, the method is 
in effect pseudo-spectral, that is, derivatives and integral operators are eval- 
uated in Fourier space, while nonlinear terms are taken in physical variables. 
In practice, the term r]%(xj,t + St) is approximated by r]j^(xj,t), and the 
equation is solved for r]i^(xj,t+St). Then this is used as a new approximation 
for the square term. This procedure is iterated until the change in the result 
falls below a certain tolerance. This scheme is second-order in time, so that 
a relatively small time step has to be used. Of course, this also means that 
the number of iterations at each time step is small, typically less than three. 

In Table 1, we record the numerical errors incurred by the time integra- 
tion of an approximate traveling wave with velocity fi = 0.8615 propagating 
for 1 and 100 periods. To find the most advantageous combinations of the 
number of Fourier modes N and the time step St, we used a computation for 
one period. We then use this combination, and integrated for 100 periods. 
The discrete L 2 -error, the difference in maximal height between the original 
wave, and the profile after 1 and 100 periods were computed. As can be seen, 
the error drops when increasing N or when decreasing h. Moreover, the fact 
that the difference in maximal height is generally smaller than the L 2 -error 
suggests that the error incurred during the time evolution is mostly due to 
a phase shift of the solution. This can also be observed in Figure [3l where 
the same traveling wave is shown after time integration for 10000 periods. 
These results also suggest that the traveling waves are orbitally stable, but 
no special investigation of this question has been carried out. 
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FIGURE 3. Solid line: approximate traveling wave <f>N for the 
Whitham equation fj, = 0.8615. Dashed line: r/jy after time inte- 
gration using (|5.4p for 10000 periods. In (a), the difference between 
<j>N and rjN is hardly visible. In this computation, N = 512 and 
St = 5e — 04. The L 2 error was 2.1e — 03, while the difference in 
height was 2.24e — 06. This and the magnification (b) suggests that 
the the error is mainly due to a phase shift. 



6. Numerical results 

In this section, we present a few numerical computations carried out us- 
ing the discretizations described above. First, we will present an in-depth 
study of some bifurcation branches of the Whitham equation, and then the 
Whitham waves will be compared to traveling-wave solutions of the KdV 
equation. 




FIGURE 4. Bifurcation branch of 27r-periodic traveling wave so- 
lutions of the Whitham equation. In (a), the entire bifurcation 
branch is indicated with dots, and compared to a curve described 
by the expressions given in (|4.7[) and (|4.6D . In (b), a close-up of 
the turning point is shown. 



The Whitham branch. In Figures Q] and El traveling wave solutions of 
the Whitham equation were shown. Figure 0] shows the bifurcation branch 
of these waves. In Figure Ufa), the whole branch is shown, and compared 
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FIGURE 5. This figure presents computations using the dynamic 
discretization with approximate Whitham traveling waves as initial 
data. The computations were done with N — 256 and time step 
St = 0.001. The initial data 4>m is represented by a dashed line, 
and the computed profile rjtq after one period is represented by a 
solid line. In (a), a wave with wavespeed c = 0.7664 is shown, In 
(b), a wave with wavespeed c = 0.7667 is shown. 



to the formula given in (|4.7p . In Figure IDJb) a close-up of the turning point 
is shown. Determining the endpoint of the branch is challenging. Here 
a combination of two tests are used. First, if a point is suspected to be 
higher than the highest branchpoint, then the solution is interpolated, and 
the Newton scheme run on a finer grid. Failure of this scheme to converge 
indicates that the point does not lie on the branch. The highest point in 
Figure Efb) corresponds to a wave for which the Newton scheme converged 
with more than 1400 gridpoints on the interval [0, 27r]. The second test 
uses the dynamic integrator. For example, Figure [5|a) shows the results of 
running the wave corresponding to the highest point on the branch shown 
in Figure 3 through the dynamic integrator for one period. Figure E{b) 
shows the next point (which is not shown in Figure |3J), but which can be 
computed with lower values of N (up to about ./V = 400) . The difference is 




FIGURE 6. A branch of approximate traveling waves for the 
Whitham equation with k = 2 is shown in (a). The profile of 
the highest wave on [0, 2ir] is shown in (b). 
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dramatic, and the time integration cannot be rescued by choosing a smaller 
time step. Nevertheless, we do not expect the numerical scheme to be able to 
approximate the wave corresponding to the very highest point on the branch. 
Indeed, comparing the maximum of the highest wave shown in Figures U] 
and [6] to the statement of Theorem 14.81 it appears that the branch of Holder 
continuous solutions should continue further than shown in those figures. On 
the other hand, our criterion for termination of the numerical bifurcation 
included failure of the time-dependent scheme to converge, and this failure 
can also be a result of instability of the computed traveling wave. Moreover, 
if the highest wave does indeed feature a cusp, then the spectral scheme will 
fail to converge, and convergence for waves close to the terminal wave will 
also be delicate. A method for computing cusped waves was proposed in pQ, 
and this might be a good approach for possilbe further studies. 

In Figure [6l a branch with k = 2 is shown. To find the highest point on 
this branch, similar methodology as for the principal branch was used. We 
note that the turning point is much farther from the end of the branch than 
in the principal branch. 

COMPARISON WITH KdV. It is of interest to see how the change in the 
linear dispersion relation affects the shape of the traveling waves. Therefore, 
attention will now be focused on the comparison of traveling waves of the 
Whitham equation to traveling-wave solutions of the KdV equation. 

First the bifurcation branches are compared. For this purpose, the wave- 
length is held fixed at 2tt. It was noted in Section [3] that the bifurcation 
points hi and /z* are n °t the same. This difference is due to the small wave- 
length 2tt. A comparison using a larger wavelength would yield a smaller 
difference in bifurcation points. Figure E[a) shows bifurcation branches for 
both KdV (solid curve) and Whitham (dotted curve) equations. As noted 
above, the Whitham branch appears to have a highest wave, after which the 
approximation no longer converges. On the other hand, the KdV branch is 
unlimited, and after a Galilean shift, this branch will approach a solitary 
wave. 

For the purposes of comparing the shape, we chose to fix the waveheight 
and wavelength of the traveling waves. Other normalizations are of course 
possib;e, such as fixing the speed or some integral measure. In Figure (7]Jb), 
(c) and (d) Whitham and KdV traveling waves with respective waveheights 
of 0.1, 0.4 and 0.485 are shown. For a better comparison, the waves are 
plotted with the minimum at zero. The small- amplitude waves are very 
similar in shape. However, for the large-amplitude waves we note that the 
Whitham waves appear narrower, with a more pronounced and sharp peak. 

7. Conclusion 

In this article, a bifurcation theory has been developed to prove existence 
of a global branch of traveling-waves solutions of the Whitham equations 
(jl.ip . Smoothness of the solutions lying on the bifurcation curve has been 
proved, and it has been shown that the solutions converge uniformly to a 
solution of Holder regularity a S (0, 1), expect possibly at the highest crest 
point (where a may be less than 1/2). The larger Holder norm blows up 
only if the highest point is attained. It was conjectured already by Whitham 
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|19] that this branch of traveling wave has a limiting wave which has a cusp. 
Unfortunately such a result has not been established, but remains an inter- 
esting open problem. A comparative local bifurcation was provided for the 
Whitham and the KdV equations, and numerical tools were used to exhibit 
properties of the principal branch (k = 1) and the secondary bifurcation 
branch (k = 2) of the Whitham equation. Finally, a comparison was made 
between the Whitham and KdV traveling waves which revealed that large- 
amplitude Whitham waves are narrower than their KdV counterparts. It 
remains to be seen which of these waves resembles the shape of a traveling 
water wave (the solution of the free-surface boundary- value problem for the 
Euler equations) more closely. 

ACKNOWLEDGMENTS. The authors would like to thank Erik Wahlen for 
pointing out an error in an earlier preprint. This research was supported 
in part by the Research Council of Norway. 
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FIGURE 7. A branch of approximate traveling waves for the 
Whitham equation with k = 1 (the principal branch) and a branch 
of traveling waves for the KdV equation are shown in (a). Com- 
parisons between the shapes of Whitham and KdV traveling waves 
are shown in (b), (c) and (d). The waves are arranged so that 
the minimum coincides with the a;- axis. The Whitham waves are 
represented by dashed-dotted curves, and the KdV waves are rep- 
resented by solid curves. The waveheight of the waves is 0.1 in (b), 
0.4 in (c) and 0.485 in (d). 
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